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SUMMARY 


When developing very high-response servosystems, the system’s linear dynamics 
can be extended by compensation and feedback techniques. This leaves nonlinearities 
and power limits as the only practical limitations of the obtainable dynamic response. 
This report demonstrates a technique for predicting the dynamic performance limits of 
a system with saturations. 

Generalized treatment for a system with a single saturatir^ nonlinearity results in 
low- and high-frequency asymptotes for upper response bounds. The bounds are veri- 
fied by comparing anal 3 dical results with frequency-response plots obtained from an ana- 
log computer model of the system. Once the amplitude dynamics are predicted with the 
limit lines, an iterative technique is employed to determine the system phase response. 

The limit -line technique is then used in conjunction with previously derived velocity 
and acceleration limits for electrohydraulic servosystems to predict the maximum dy- 
namic performance obtainable for an electrohydraulic servosystem containing a single- 
stage servovalve. The servovalve used in this system was a single-stage electrohydrau- 
lic design that provides linear dynamic capabilities beyond the range of conventional two- 
stage electrohydraulic servovalves. By deriving relationships that define the upper 
bound on system performance, those system parameters that affect these bounds can be 
identified. The resulting design criteria can be used to optimize the system design. 

The predicted maximum performance for the servosystem is compared with experi- 
mentally obtained response data. The results demonstrate that the limit-line predictions 
agree quite well with the data. The optimized system exhibited an amplitude response 
that was flat to within ±3 decibels of the commanded amplitude for frequencies to 530 
hertz for a peak-to-peak displacement of 0. 127 centimeter and a load mass of 1. 1 kilo- 
grams. 


INTRODUCTION 

Electrohydraulic, servoactuated air valves have been used extensively as research 
control and perturbation devices on turbojet engines and supersonic inlets (refs. 1 to 3). 
As perturbation devices, they provide a means for creating pressure or flow perturba- 
tions in the normal flow processes of engine components. Dynamic analysis of the re- 
sulting signals provides transfer-function models that can be used to describe the sys- 
tem (refs. 4 and 5). When these valves are placed in control loops, they can be used for 


stabilizing unstable flow phenomena or for stud 3 ^g control techniques (ref. 6). To gen- 
erate wide -bandwidth dynamic data or to implement fast control loops, these valves 
must be capable of responding to high-frequency dynamic signals. 

In designing such a servosystem, the designer usually considers the linear dynam- 
ics of the electrohydraulic servosystem and finds ways to optimize the linear response 
through available control-system techniques (ref. 7). However, once the linear dynam- 
ics have been extended through compensation, the servosystem becomes limited in its 
dynamic performance by such nonlinearities as signal saturations or limited available 
hydraulic power. A previous analysis (refs. 8 and 9) of servosystems containii^ a two- 
stage electrohydraulic servovalve was performed to determine the effect that limitii^ 
the hydraulic power available to the actuator would have. That analysis derived acceler- 
ation and flow limits on the dynamic response of the hydraulic actuator piston for peak- 
power-transfer operation. Such performance limits are nonlinear and tend to degrade 
the servoactuator’s linear dynamic performance for large -amplitude motions. The rela- 
tionships that define the limits relate the amplitude and frequency of the output sinusoid. 
On a log-magnitude -versus-log -frequency plot, these limits are straight lines that de- 
fine the upper performance bound. This technique has been successfully used to design 
overboard bypass doors for supersonic inlets (refs. 2 and 3). It has also been used to 
design a high-response airflow valve for generating upstream airflow perturbations and 
distortions in a turbojet-engine altitude test chamber (ref. 1). 

This report extends the anal 3 d:ical technique derived in reference 9 to the case where 
a saturating nonlinearity exists in the generation of a position error signal. The system 
described herein was designed with a single-stage electrohydraulic servovalve consist- 
ing of a four-way hydraulic spool valve centered with springs and driven by a high- 
current coil similar to an audio speaker driver. The sii^le-stage servovalve was in- 
tended to provide better linear dynamics than a conventional two -stage servovalve. By 
using feedback compensation around the spring-mass spool system, the linear frequency 
response can be extended to the point where system nonlinearities become the major 
determinant of servosystem dynamic response. In this system, saturation occurs at the 
output of the power amplifier that drives the coil. Therefore, a technique is needed to 
readily predict how this saturation affects the system response. Such a technique would 
also provide a means for selecting such design parameters as spring rate and actuator 
piston area. 

Althou^ the analysis reported herein is used for electrohydraulic servosystems, 
its application to other servosystems containing saturation is apparent. To provide a 
more useful approach, this report derives a simple, generalized technique for deter- 
mining how saturating nonlinearities affect the maximum obtainable frequency response 
for a typical feedback control system. The analysis derives two lines on the log- 
amplitude -versus -log-frequency plot that define system operation in the nonlinear re- 
gion. The first line is determined by assuming that the saturating device puts out a sine 
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wave with a peak amplitude equal to the saturation level and then applying this amplitude 
to the low- and high-frequency asymptotic approximations of the system dynamics. This 
line, called the linear-nonlinear transition bound, defines a locus of points where the 
system initially crosses from total linear operation to operation with saturation. The 
second line is derived in a similar manner, except that the output of the saturating de- 
vice is assumed to be a square wave and the amplitude of the square wave’s fundamental 
frequency is applied to the system dynamics. This line, called the maximum - 
performance limit line, provides a locus of points beyond which system operation is the- 
oretically not possible. 

These analyses use asymptotic approximations to system dynamics; therefore, the 
locations of the bounds defined by these lines may be in error near breakpoints in the 
asymptotes or at system resonant frequencies. For most systems with well -damped 
poles, these simple lines are adequate for performance predictions and, as such, pro- 
vide a very quick and easy approach. However, if the system being analyzed is known 
to have lightly damped poles, the limit lines must be derived by considering the exact 
magnitude of the system d3mamics at all frequencies. This results in a more complex 
relationship and a limit -line plot that is not linear on the log-amplitude -versus -log- 
frequency plot. An example of such a system is provided. 

The limit lines (both simple and complex) for the example system are combined with 
the system linear dynamics in an iterative analysis to determine the system phase -angle 
response. This results in a family of phase curves with the input command amplitude as 
the parameter. 

The analysis of the example system is presented and compared with an analog com- 
puter model. The limit-line technique is then applied to the design of an electrohydrau- 
lic servosystem driven by a single-stage servovalve. The saturation technique is used 
to predict how coil voltage and current saturations affect system dynamic response. 
These are combined with the actuator velocity and acceleration limits to provide a com- 
plete prediction of system performance. The actual servosystem was tested by using 
frequency sweeps at various commanded amplitudes. The test results are presented in 
the form of nonnormalized amplitude and phase plots and are compared with the pre- 
dicted performance. 

This work supported the design of a specific actuator. The actuator was used both 
as a key control element in a stabilizing control system and as a high-frequency pertur- 
bation device for dynamic studies. However, the fundamental assumptions and relation- 
ships derived should be applicable to any servosystem containing similar nonlinearities. 
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DERIVATION OF SATURATION PERFORMANCE LIMITS 


Simple Limit Lines for a System with Saturation 


To provide a somewhat generalized treatment of a servosystem operating with its 
error signal in saturation, we analyzed the simplified unity -gain feedback system shown 
in figure 1. All symbols are defined in appendix A. The error signal E was generated 
by differencing the command signal R and the output signal C. This error was multi- 
plied by the gain and limited by the saturation block (labeled to provide a force 
signal F. The system dynamics are represented as a second-order lag with a free inte- 
grator and a steady-state gain Kp. 

To simplify the analysis of the system when it is operatii^ in saturation, we first 
assumed that only the low-frequency dynamics are important (| as^ + bs| < 1) and then 
that only the high-frequency dynamics are important (|as^| > jbs + l|). This resulted 
in solutions that asymptotically approximate the d 3 uiamics of the second-order pole. 

Low-frequency linear performance bound . - From figure 1, the transfer function 
that relates the output to the force for low-frequency dynamics can be written as 


C(s) _ ^p 
F(s) s 


( 1 ) 


If the force signal is a sinusoid whose peak amplitude just touches the saturation level (or 
I F(s)| = letting s = jw and taking magnitudes give 

( 2 ) 

(V 


If this relationship is plotted as | C | versus on a log-log plot, the resulting straight 
line is the dynamic bound between linear and nonlinear operation. Figure 2 presents a 
plot of this line (the low-frequency segment of the dashed lines) for three values of Kp. 
The numerical parameters chosen for this system are presented in table 1. 

High-frequency linear performance bound . - The output transfer function for high- 
frequency dynamics ^ig. 1) is 


C(s) _ S 


P(s) 


as 


(3) 


Again, letting |f(s)| = Fg^t s = jw give 
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( 4 ) 


K F 

|c| = ... P sat . 

which is the Bode -plot locus of points for the high-frequency transition from linear to 
nonlinear operation. These lines are the high-frequency segments of the dashed lines in 
figure 2. They have a slope of -3 decades per decade; the low-frequency line has a 
slope of -1 decade per decade. 

Limit lines . - Although the previously derived lines are convenient for defining the 
transition from linear to nonlinear system performance, they do not provide information 
as to the maximum achievable dsmamic performance with this system. Maximum per- 
formance is obtained by assuming that the error signal E is at a sufficiently large am- 
plitude to make the saturating -device output F approach a square wave with an ampli- 
tude that is set by the saturation level. The fundamental frequency of the square wave is 
the same as the excitation frequency. 

Describing-function theory (refs. 10 and 11) defines the output from a saturating de- 
vice as 


f(t) = Bj sin ojt + ^ 

n=3, 5, 


B^ sin no)t 


for an input of 


Kg e(t) = Kg I E I sin cot 


(5) 


( 6 ) 


where 


B 



(Ke|E|) 



^sat 

KelEl 



(7) 


To approximate this signal, we considered only the fundamental frequency of the force 
signal f(t). The describing function coefficient B^ was determined by taking the limit 
of equation (7) as the error-signal amplitude | E | approached infinity. 


B 



( 8 ) 


This gives 




( 9 ) 
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Using this input for deriving the plant high- and low-frequency ou^t amplitudes results 
in 

Low-frequency limit: 


4K F__4. 
C = P 


ffO) 


( 10 ) 


Hi^-frequency limit: 


|C| = ■ g (11) 

Q 

Traco'^ 

The lines defined by equations (10) and (11) represent the maximum output amplitude 
obtainable for a given frequency. Figure 2 shows the maximum -performance lines 
(solid) for three values of Kp. Essentially, they are the linear -nonlinear transition 
bounds shifted in magnitude by a factor of 4 /tt. The region between the dashed and solid 
lines defines the operatii^ region where the saturating -device output is a clipped sine 
wave - the nonlinear performance region. In all cases, the low-frequency lines inter- 
sect the high-frequency lines at the natural frequency of the second -order pole (100 Hz). 
The solid lines are referred to as the system -performance limit lines throughout the re- 
mainder of the report. 

Analog computer model of example system . - To verify the derived limit lines, we 
implemented the system of figure 1 on an analog computer and subjected it to frequency 
sweeps with various input -commanded amplitudes | R [ . The tests were performed with 
the second-order -pole damping ratio C at 0. 7 for the three values of and Kp 
shown in table I. The product K K was kept constant to maintain the same closed- 

e P 

loop linear dynamics for all cases. 

The results of these tests are shown in figure 3 with the respective limit lines drawn 
in. System ou^ut amplitude | C | is plotted against frequency on a log-log plot. Figure 
3(a) (K = 10, K = 50) shows that the large -amplitude response (|r| =5) was degraded 
and followed the limit line for frequencies above 16 hertz. However, when the forward- 
loop gain K K was redistributed to place more gain downstream of the saturating de- 
vice, the limit lines moved up in amplitude and were less of a restriction on large - 
amplitude performance, as shown in figure 3(b) for gains of K = 5 and K = 100. 

e p 

When Kg was decreased to 1 and Kp was increased to 500, the system operated lin- 
early for amplitudes |r| of 0. 5, 1.0, and 2.0, as shown in figure 3(c). For all cases 
shown here, the actual system response was somewhat lower than the limit lines at the 
100-hertz natural frequency of the open-loop pole. This prediction error is due to the 
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asymptotic approximations to the second-order-pole dynamics that were used to derive 
the limit lines. 


Derivation of Complex Limits 

The previous derivation used asymptotic approximations to the second-order pole in 
the plant dynamics. Often a servosystem may contain a lightly damped pole, in which 
case asymptotes display considerable errors when used to approximate the system dy- 
namics. To analyze such a system, the limit lines must be derived by considering the 
exact magnitude of the system response at all frequencies. 

Analysis of saturation followed by a lightly damped pole . - For the system in fig- 
ure 1, the transfer function for the dynamics between C and F can be written as 

( 12 ) 

s(as2 + bs+« 

If sinusoidal excitations are assumed, s = j can be substituted and the saturating- 
device output can have a peak amplitude equal to the saturation level (| F(s)| = 

Then the output amplitude is 


Id Vsat (j3) 

W ^ a^«4 + (b^ - 2a)w^ + 1 

This relationship defines the linear -nonlinear transition bound, as before, but in a more 
complex form. Similarly, the maximum -performance limit line is found by assuming 
the saturating -device output to be a square wave and taking the fundamental-frequency 
amplitude to be 4 /tt times the saturation level. Therefore, the complex limit line is 

|C| = — (14) 

ttO) y^a^w^ + (b^ - 2a)<^ + 1 

The complex limit line is plotted in figure 4 for various damping ratios and is compared 
with the asymptotic limit line for the case where = 5 and Kp = 100. From this fig- 
ure, we can conclude that the asymptotes give good results for systems with damping 
ratios near critical damping but underestimate the response for lightly damped systems 
and overestimate it for highly damped systems. 

Analog computer model of lightly damped example system . - Again, to verify this 
analysis, we ran a lightly damped case on the analog model. However, when trying to 
place a lightly damped pole in the forward path of a closed-loop system, it is difficult to 



close the loop and obtain good relative stability. In a practical system, an attempt is 
usually made to provide additional damping by closing a velocity feedback loop around 
the lightly damped pole. The electrohydraulic servosystem to be discussed later is such 
a system. To provide an analysis relevant to that case, and also to allow for compari- 
son between the ? = 0. 2 and ^=0.7 results with the example system, we included a 
velocity feedback loop, as shown in figure 5. The damping ratio of the second-order 
pole was set to a dampii^ ratio of 0.2 (P3), and the velocity feedback loop P5 was 
placed outside the saturation device to create an effective damping ratio of 0. 7 when the 
system operated linearly. Then, when the system went into hard saturation, the loop 
would be effectively opened and the damping ratio would drop to approach 0.2. 

The result of this test is presented in figure 6, along with the complex limit line and 
the linear-nonlinear transition bound. The large -amplitude responses for j r| of 5 and 
9 are relatively flat until they intersect the limit line, after which the response curves 
follow the limit line closely. The sweep that was done at an amplitude |r| of 2 does not 
reach the limit line at all. But it does cross the linear -nonlinear transition boimd, and 
the partial saturation distorts the response somewhat from the linear response. In 
terms of linear dynamics, the system presented in figure 6 is essentially the same as 
the one presented in figure 3(b). But comparing these two figures shows that the system 
that has a lightly damped pole following the saturation can extend its large -amplitude re- 
sponse beyond that obtained with a highly damped system. For example, the system re- 
sponse for an |r| of 5 in figure 3(b) is down 3 decibels at a frequency of 43 hertz, but 
for the response in figure 6 the -3-decibel point is at 103 hertz. This implies that, in 
designing such a system, lightly damped poles should be placed downstream of the sat- 
urating device and that the needed damping should be provided with a velocity feedback 
loop around the saturating device. 


System Phase Response 

Often the system to be designed will be used as an element in a larger system. In 
this case, a phase-response prediction is necessary to determine overall system stabil- 
ity. Saturation creates additional phase beyond that found in linear systems, and the 
added phase is a function of the system command amplitude. Solving for the phase is 
difficult since either a piecewise linear model in the time domain or a detailed analog 
computer model is needed. 

Phase determination from limit lines. - This section briefly describes a technique 
for determining phase in a saturating system. The technique uses the limit-line equa- 
tions and the system linear dynamics. A more detaited discussion of the technique and 
its digital implementation and the digital program listings are presented in appendix B. 
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Basically, the technique assumes that the system operates linearly with the satura- 
tion characteristic replaced by a gain term, This gain is set to unity for all oper- 

atii^ conditions where the system linear output amplitude is less than any one of the 
limit-line amplitudes at a given frequency. When the linear solution for amplitude, at a 
given frequency and command amplitude, exceeds any of the limit-line amplitudes, the 
value of is decremented iteratively until the system output amplitude equals the 
amplitude of the lowest limit line. The resulting linear system with a decreased 
forward-loop gain is assumed to be a valid model for the system at that particular com- 
mand amplitude and frequency. The resulting phase between command signal and output 
signal is assumed to be the correct system phase at that point. 

Analog computer model comparison . - The results for this scheme and the analog 
model data are presented in figure 7(a) for a damping ratio of ^ of 0. 7 and a plant gain 
Kp of 100. The solid lines represent the analog-model phase response and the data 
points are the results from the iterative phase technique. Althoi^h the simple limit- 
line equations used to determine these points were asymptotic approximations, the phase 
results agree well with the analog response. The additional phase contributions caused 
by larger amplitude excitations can clearly be seen. 

The case for the system with an open-loop damping ratio of 0. 2 is shown in figure 
7(b). The iterative phase technique in this case used the complex limit-line equation and 
took into accoimt the effective opening of the velocity feedback loop when the device goes 
into saturation. Additional detail is provided in appendix B. Again good agreement was 
obtained between the results of the iterative technique and the analog model. 


Multiple Saturations 

As the complexity of a control system increases, the possibility of having several 
amplifiers or devices that can enter into saturation becomes more likely. K one satura- 
tion is predominant, the others can be ignored and the previous analysis can be used. 

But many times, dynamic elements are present between the saturating devices. Which 
device enters into saturation first is a function of amplitude and frequency. The previ- 
ous analysis can be used to determine which device is limitii^ the system performance 
at any given operating condition by ignoring all but one of the saturations and solving for 
its limit lines. The same is done for each saturation in succession. This generates a 
set of limit lines that can be compared in order to determine the most restrictive limit. 
As frequency increases, the system may shift from one limit line to another, as in the 
case of the low- and high-frequency limit lines in figure 2. In all cases, thi line that 
most severely limits system response at a given frequency identifies the device that is 
in saturation and thus is limiting performance. 
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No rigorous proof is provided for this superposition theory of limit lines. But, as 
shown in later sections of this report, this technique does agree well with the escperi- 
mental results. 


APPLICATION OF LIMIT -LINE TECHNIQUE TO DESIGNING A HIOI- 
RESPONSE ELECTROHYDRAULIC SERVOSYSTEM 

The limit-line technique for saturation was developed because of problems in de- 
signii^ servosystems containing a single-stage electrohydraulic servovalve. This sec- 
tion of the report describes a high-response servosystem, derives a set of performance 
limit lines based on saturation analysis and power limitations, and compares the e^)eri- 
mental results with the predicted d 3 mamic response. The system has both a lightly 
damped resonance and multiple saturations that must be treated in the analysis. 


Description of Airflow-Valve Servosystem 

An airflow valve was des^ned to provide a variable and controllable orifice area for 
studying the dynamics and controlling the stability of an axial compressor. A schematic 
diagram of the actuator servosystem used to control the valve area is shown in figure 8. 
Since the discussion of this report is limited to the servosystem, the slotted cylinder 
that forms the airflow-valve movii^ element is treated merely as a load mass on the 
actuator piston. 

Experimentally, the actuator piston position was measured with a linear variable 
differential transformer (LVDT). Its output signal was fed into an electronic controller, 
where it was differenced with a precompensated command signal and a set-point signal. 
The precompensator, a second -order lead-lag circuit, was used to extend the linear re- 
sponse of the servosystem. The error signal was amplified and applied to a power am- 
plifier that drove the coil on the electrohydraulic servovalve. An Internal feedback loop 
provided current feedback to help extend the amplifier-coil dynamics beyond the 
resisitive -inductive dynamics of the coil. 

The electrohydraulic servovalve is a single-stage, four-way spool valve positioned 
mechanically by the sum of the forces exerted by the springs at either end of the spool 
and the current-driven coil. The spool valve opens a port to hydraulic supply and return 
pressures intS the actuator piston chambers in order to position the piston. The second- 
order dynamics of the spool mass suspended by springs is very lightly damped and re- 
quires velocity feedback to remain stable. This feedback is provided by the linear veloc- 
ity transducer (LVT) mounted at the left end of the servovalve spool in figure 8. The 
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output signal of the velocity transducer is fed into a feedback compensator, which pro- 
vides essentially an integral-plus -proportional feedback, to provide added damping and 
to increase the sprii^-mass resonant frequency. 

A detailed block diagram of the servosystem is provided in figure 9. The system is 
divided into major elements identified by the dashed lines. The controller is an analog 
circuit that generates the signals to be fed to the power amplifier. The second-order 
precompensator is shown actir^ on the command signal. A slew -rate limiter is used to 
prevent the compensated command signal from requestii^ very high velocities at large - 
amplitude excitations. The slew -rate -limiter output is summed with the actuator posi- 
tion feedback and the set point. The error is multiplied by K and fed to the power 
amplifier. The velocity signal from the servovalve spool is modified by a band- 

pass filter to provide a high -pass -filtered and integrated spool velocity signal. This 
signal allows the spool resonant frequency to be increased without decreasing the steady- 
state, effective forward-loop gain. 

The power amplifier is modeled as a high-gain amplifier whose output Op^ 

can saturate. The power -amplifier voltage output is applied to the first-order lag of the 
servovalve coil formed by R and L to generate a coil current i that also has a 
saturation level. This saturation is built into the amplifier as an adjustable limit whose 
maximum is determined by the coil thermal operating limits. The current is measured 
by the voltage drop across the series resistor Rj and is fed back into the power ampli- 
fier. Since the current-measuring resistor is in series with the coil, the voltage drop 
across it is subtracted from the power-amplifier output voltage before it is applied to the 
coil impedance. Coil current passing through the permanent-magnet field creates 
an excitation force f that is differenced with the spring force K x to accelerate the 

servovalve spool mass M . The spring-mass spool mechanical system is modeled as a 

s 

pure second-order system whose output is spool-valve position x . 

s 

The spool position is the input to the servovalve orifice equations that determine the 
hydraulic flows and pressures in the actuator cavities (refs. 12 and 13). The spool ori- 
fice areas are computed by assuming symmetry and applying an open-center bias area 
Aq. Four orifice equations are used to solve for flows that are applied to the actuator 
volume dynamics in order to determine the chamber pressures Pj,P 2 - The pressure 
difference p^^ across the piston accelerates the piston load mass M^, and a double in- 
tegration provides piston or airflow-valve position x^. 

To analyze the dynamic performance limits of this system, we must consider the 
saturations of the coil voltage e^^ and coil current i^, the lightly damped resonances 
formed by the servovalve spool dynamics (x to f„), and the actuator hydraulic volume 
dynamics (x^ to x^). The physical parameters that were used to model the airflow- 
valve servosystem are listed in table n. 

The airflow valve mounted on a test stand is shown in figure 10. The servovalve 
was rotated 90° relative to its position in figure 8 to prevent case vibrations caused by 
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piston motion from coupling into the spool position. This could create an instability. 

The servovalve was attached directly to the actuator block to minimize the hydraulic vol- 
ume of the actuator control ports. 


Actuator-Piston Acceleration Limit 

The limit lines derived in this section and the next are similar to those derived in 
reference 9 and pertain to the power-drive limitations imposed by the finite supply pres- 
sure and flow orifice areas of the servovalve. The load that the actuator piston drives 
against is assumed to be a pure inertia load set by the mass of the slotted, cylin- 
drical airflow valve and the actuator piston. The force balance on the piston does in- 
clude a damping term (fig. 9) that is very small relative to the pressure forces and is 
neglected for this analysis. The resulting equation is 

= ^a^Pl ■ = Va 

The dynamic performance limit for the actuator piston is obtained by assumir^ that the 
pressure drop across the piston is at its peak value. Under ideal circumstances, the 
maximum force available to accelerate the actuator piston is estimated by assuming one 
chamber to be at full supply pressure P_„ and the other chamber to be at zero pres- 
sure. However, a better approximation to peak pressure is obtained by assuming the 
piston pressure -flow characteristics to be operating on an ellipse that delivers peak 
power to the actuator. A complete derivation of this curve and the resulting peak pres- 
sures and flows are presented in appendix C. The resulting peak pressure is 

Pa. peak = 

Assuming sinusoidal output motion (x^ = |X^| sin cut) and equating peak values for equa- 
tion (15) give 

, , « 2\/2A„P_„ 

X cu2 = (17) 

^ 3M^ 

which is the expression for the actuator acceleration limit. 
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Actuator Piston Velocity Limit 


The actuator piston velocity is essentially proportional to the hydraulic flow enter- 
ing the piston chambers. If the fluid is incompressible, 

qg = qg = % (18) 

Based on this assumption, the flow balance for one volume of the actuator piston cham- 
ber would be 


■is = \\ 


(19) 


The maximum possible velocity is determined from the maximum servovalve flow. To 
find what the maximum flow would be, we use the peak -power elliptical load line of ap- 
pendix C. That is, 

(q«) = V2 Qr (20) 

® peak ^ 


where is the servovalve rated flow at the peak -power transfer point and is defined 
as the flow through the servovalve orifices when the load pressure on the piston is 


P 


a 



( 21 ) 


and the spool valve is at its maximum open condition. The flow through the servovalve 
orifices for a positive displacement of x can be written as 






^^su “ 


( 22 ) 


= ’^'^s^s^d 



( 23 ) 


, From equation (18) we can set equation (22) equal to equation (23) to get 

^su = Pi + P2 (24) 

Using the assumption of matched symmetrical orifices, we derive the following relation- 
ship for flow: 

‘Is = ’^Vs^d |/^(Psu-Pa> (25) 
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(ref. 12, pp. 81-83). Substituting equation (21) into equation (25) gives the rated flow 


Qr^^^^sl^sl ^d 

max 


su 

3p 


( 26 ) 


where |x | is the maximum spool displacement amplitude. 

“ max 

Substituting equations (20) and (26) into equation (19) and considering only amplitudes 

give 





(27) 


where 



(28) 


This equation is the velocity limit for the airflow-valve servosystem. However, the 
maximum displacement of the servovalve spool is not a constant value. It depends on the 
saturation levels of the current and voltage signals that drive the spool coil. The term 
can be viewed as the linearized gain between the spool position and the serovalve 
flow. 


Coil-Current Limit 

From the block diagram of figure 9, we can write the transfer function for servo- 
valve spool position to coil current as 


For only low frequencies. 


Xs(s) 



^s 2 °s , 
s + s + 1 


K 


sp 


K 


sp 




(29) 


( 30 ) 
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Assuming a sinusoidal coil current whose peak amplitude is at the saturation level and 
considering only magnitudes give 


X I = 


Ksp 


( 31 ) 


which is the spool displacement amplitude at which the system crosses the linear - 
nonlinear performance bound. As was shown earlier, the maximum possible spool dis- 
placement at low frequencies is given by 


c c, sat 


max 




(32) 


sp 


Substituting equation (32) into equation (27) gives the low-frequency coil-current limit. 

4K,K I, 




c c, sat 


(33) 


Returning to equation (29) and considering only high frequencies yield 

Xg(s) 


Ie(s) 


MgS‘ 


(34) 


Substituting s = jw, taking magnitudes, and letting |lj.(s)| = 1^ give 


I I 2 ^c 


M 


c, sat 


(35) 


s 


for current peaks at the saturation level. Solving for |Xg|, applying the 4 /tt factor, 
and substituting into equation (27) give the high-frequency coil-current limit. 

1^.1 


0,3 = 


4KzKeIc,sat 


77M„A„ 

s a 


(36) 


Coil -Voltage Limit 


The presence of coil-current saturation is ignored in this section, and the effect of 
coil -voltage saturation is investigated. The linear transfer function between the spool 
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position and the coil voltage includes the dynamics of the servovalve coil and the sprir^- 
mass dynamics of the spool. Since this analysis treats the case of saturation, the spool 
velocity feedback loop is assumed to be effectively open. The transfer function is 


V“> 


K. 


K „(R, + R„) 
sp' 1 c' 




s + 1 


K 


sp 


K 


sp 


(Rf + Rj. 


s + 1 


Assuming only low frequencies gives 


Epa(s) 


(37) 


(38) 


If we use the same techniques as in the previous section, the low-frequency limit line 
becomes 


|x I ... - sat 

'' ’'V''! * ®c)Aa 

Similarly, the high-frequency transfer function is 

Xs(°) _ Kq 

and the high-frequency coil-voltage limit line is 

|X I r.A _ '^^Z^c^pa, sat 


(39) 


(40) 


(41) 


The limit -line equations for the airflow-valve servosystem were determined from 
the parameters of table n and are presented in table m in summary form. Figure 11 
provides a plot of these equations, showing the linear -nonlinear transition bounds as 
dashed lines and the maximum -performance limit lines as solid lines. The coil-voltage 
limits are beyond the coil-current limits for almost all frequencies and, therefore, they 
are only a factor for frequencies beyond 900 hertz. The acceleration limit line is lo- 
cated just beyond the coil-current limit lines. 
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Using Limit Lines to Design the System 


Inspection of the limit-line equations of table HI reveals that a common parameter 
is actuator piston area A_. By selecting this area, the position of tiie limit lines can be 
selected. The magnitude of the acceleration limit line is a direct function of piston area, 
but the coil current and voltage limits vary inversely with piston area. The movement 
of a t3Toical set of coU-current and acceleration limit lines for decreasing piston area is 
shown in figure 12. The coil -current lines move up in amplitude and frequency while the 
acceleration limit line moves down. To select a piston area that would tend to maximize 
the system response, we equate the point where the low- and high-frequency current- 
limit lines intersect with the acceleration limit line. 

The intersection of low- and high-frequency current limits is found from 


Kl = 


4K,K I 4K,K„r 

L C C y S3.W __ L C Cy S2tt 


sp a 


(42) 


TrMsAa^^" 


which results in 


CO 


♦ 



(43) 


which, of course, is the resonant frequency of the spring-mass spool dynamics. Eval- 
uating the acceleration limit at the frequency defined in equation (43) gives 


or solving for A* 

a 


X.I = 


4K?K^Ic,sat _ 


_ 2 |/^ 




3/2 





(44) 


(45) 


which is the actuator area that will make the limit lines intersect, as shown by the 
dashed lines in figure 12. From the parameters in table n, the optimum area is 2.96 
square centimeters. The actual area used, 3.20 square centimeters, was a design com- 
promise based on piston-seal sizing considerations. We decided to use a larger than 
optimum area because this made the coil-current limit the predominant limit. Since the 
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dynamics downstream of the current saturation contained a lightly damped resonance, 
the actual system response will probably exceed the coil-current limit line near the res- 
onance. 


Complex Limit Lines for Airflow -Valve Servosystem 


As was discussed in earlier sections, the presence of lightly damped poles down- 
stream of the saturating device requires a more exact solution for the limit lines, one 
that uses transfer-function magnitudes rather than asymptotes. Figure 9 shows that, 
besides the resonance of the servovalve spool spring-mass, there exists a second-order 
resonance formed by the stiffness of the hydraulic fluid /3 driving the piston mass M^. 
This resonance was ignored when the limit lines were derived because the actuator vol- 
umes and the piston and airflow -valve masses were sized to set this resonant frequency 
beyond the range of interest. However, practical design considerations lowered the res- 
onant frequency to 472 hertz, which is near the 3 15 -hertz resonance of the servovalve 
spool. Therefore, when the transfer function used to derive the complex limit lines is 
being written, the piston hydraulic resonance should also be considered, as shown by the 
simplified block diagram of figure 13. 

The following derivation will be carried out for the coil-current limit since it is the 
most restrictive of the saturating limits. The coil-voltage limit could be determined 
similarly by including the coil impedance between and i^. 

The original derivation of the actuator velocity Umit basically used the same block 
diagram with the transfer function of x^ to qg replaced by 1/A^s. The transfer func- 
tion from Ig(s) to Ky(s) can be written as 


Xy(s) 

Ic(s) 


SE_A. 


'^s 2 ^s 
sf — ^s ^+ — — s + 1 


‘sp 


K 


sp 


2^ 

s + 


5m 


M.. 


(46) 




s + 1 


The spool flow to pressure gain K^p q was determined from the relationship in refer- 
ence 12 (p. 97). The denominator terms can be simplified by rewriting them in natural 
frequency and dampii^ ratio form 
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(47) 


where 


\(s)_ 


^sp^a 


s/li.s2 + !fis+lV-i-s2 + !ills+r 


tO)‘ 


0) 


ns 


ns 




nh 



(0 


nh 


D, 


= 




‘"nh = 2A 


3p,o-i/^ 
2 Aa r 


(48) 


(49) 


(50) 


(51) 


Assuming the input i^ to be the fundamental component of a square wave and the output 
to be a sine wave with a zero-to-peak amplitude of |X^| and taking the magnitude of 
the transfer function by replacing s with jw give 



4K,K I 


I c c, sat 
^ ^spAa 



+ j2?. 


O) 


iW 


ns> 



+ j2? 



(52) 


which is the equation for the complex coil-current limit line. The limit line is plotted in 
the following section for the airflow-valve servosystem. 


RESULTS AND DISCUSSION 

Experimental data for the dynamic response of the airflow-valve servosystem were 
obtained by sweep frequency -response techniques at various excitation amplitudes. The 
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frequency responses were plotted with the log of nonnormalized amplitude on the vertical 
axis and the log of frequency on the horizontal axis. The results are shown in figure 14 
for excitation amplitudes of 5, 10, 20, 40, 80, and 100 percent of full stroke (full 
stroke = 0. 318 cm, zero to peak). For comparison, the acceleration limit line and the 
complex current-saturation limit line are plotted on this figure. The amplitude response 
of the servosystem moving a 1. 1-kilogram mass was flat to within 3 decibels of the low- 
frequency amplitude to 530 hertz for a 20 -percent amplitude signal. This response de- 
graded to 200 hertz for an 80 -percent amplitude signal. 

The output was limited at large amplitudes by the slew -rate limiter used to precom- 
pensate the command signal in order to avoid excessive velocities. At frequencies be- 
tween 300 and 550 hertz, the system resonated through the current limit asymptotes to 
the acceleration limit line. For frequencies beyond 550 hertz the amplitude dropped off 
along the complex current limit line. This amplitude result agrees well with the pre- 
dicted performance derived from the limit-line analysis. The small-amplitude response 
at 10 -percent amplitude is essentially the linear response of the system for frequencies 
to about 600 hertz. At 5 -percent amplitude, stiction tends to decrease the resonant peak. 

Figure 15 presents the phase response of the airflow-valve servosystem for 20- and 
100-percent excitation amplitudes. The system crossed the -90° point at 270 and 140 
hertz for 20- and 100-percent amplitudes, respectively. The data points shown on the 
curves are the phase prediction from the iterative analysis that was based on the pres- 
ence of the acceleration limit and the complex current limit line. The slew -rate limiter 
was also modeled to achieve agreement with the 100 -percent-amplitude result. 


CONCLUSIONS 

An analysis for determining the dynamic response limitations for a servosystem 
containing saturations has been presented for a simple second-order system and for a 
practical electrohydraulic servosystem design. The boundary between linear and non- 
linear performance is determined for saturation by assuming a peak amplitude out of the 
saturating device equal to the saturation level and by applying this amplitude to the linear 
dynamics between the saturating device and the system output. The maximum- 
performance limit line is the same as the nonlinear bound line, but with the system out- 
put amplitude multiplied by 4/7t. 

Having applied the limit-line criteria to these systems, we can draw the following 
conclusions: 

1. When designing a system with saturation, it is desirable to distribute as much of 
the forward -loop gain as possible downstream of the saturating device. 

2. The use of asymptotes to approximate linear performance, is reasonable for sys- 
tems with damping ratios of 0. 5 to 1.0. However, systems with lightly damped, open- 
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loop poles downstream of the saturation may require a more detailed analysis near the 
resonant frequency. This analysis shows that the system performance can exceed the 
asymptotic maximum -performance limit line because of the system's high gain at the 
resonant frequency. 

3. Decreasing the damping on the open-loop pole and creating the needed damping by 
using a feedback loop around the limiting device can extend the response of a system 
near the open-loop resonant frequency. This technique results in a response somewhere 
between the asymptotic limit lines and the detailed complex limit line. 

4. For systems with more than one saturation, the limit lines can be obtained by 
treating one saturation at a time and assuming the other saturating devices to be in their 
linear operation regions. The most restrictive limits are then used to define the 
maximum-performance region. 

5. The optimum actuator piston area for an electrohydraulic servosystem can be 
determined from the limit-line equations to provide the maximum system performance. 

6. If the design specifications for a servosystem stipulate maximum allowable ve- 
locities or accelerations, parameters can be chosen to position the limit lines to assure 
that these specifications are not exceeded. 

The iterative technique of combining the predominant limit lines with the system 
linear response to determine system output phase response provides good agreement 
with actual system response. The assumption that operating in saturation opens any 
feedback loops around the saturating device is an oversimplification that may require 
further investigation. Since the linear feedbacks are functional over the portion of the 
waveform period when zero-error crossover occurs, a modified set of relationships 
could be developed to treat this case. Although the scope of this work precluded enter- 
ing into this type of detail, it is suggested for future studies. 

The fast-response electrohydraulic servosystem designed with the limit -line tech- 
nique provided good dynamic response. It moved a mass of 1. 1 kilograms through a 
peak-to-peak displacement of 0. 127 centimeter at 530 hertz (±3 dB), with a phase angle 
of -90° at 270 hertz. This response agreed well with the predicted results from the 
limit-line analysis. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, March 30, 1979, 

505-05. 
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APPENDIX A 


SYMBOLS 

A area, cm^ 

2 

a example -system transfer -function coefficient, sec 

B describing -function coefficient (pseudo force signal) 

b example -system transfer -function coefficient, sec 

C orifice discharge coefficient 

C(s) example -system output signal (frequency domain) 

c example -system output signal (time domain) 

D damping factor, N • sec/cm 

d diameter, cm 

E(s) example -system error signal (frequency domain) 

e airflow -valve servosystem voltage, V 

e(t) example -system error signal (time domain) 

F force level, N 

F(s) example -system pseudo force signal (frequency domain) 
f(t) force (time domain), N 

I(s) current (frequency domain), A 

i current (time domain), A 

K gain 

servovalve coil-force-to-coil-current gain, N/A 

Kjp Q servovalve null flow -pressure coefficient, cm®/N • sec 

servovalve linearized flow to maximum spool position for peak-power transfer 
point, cm^/sec 

power amplifier gain, V/A 
pa 

spring rate, N/cm 
sp 

K proportional gain for spool velocity feedback, V/V 

ps 

K. integral gain for spool velocity feedback, sec 

IS 

spool velocity transducer gain, V* sec/cm 
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airflow-valve position transducer gain, V/cm 

L inductance, H 

p 

M mass, N* sec /cm 

p 

P pressure (steady state), N/cm 

p 

p pressure (instantaneous), N/cm 

O 

Q volumetric flow (steady state), cm /sec 

O 

q volumetric flow (instantaneous), cm /sec 

R resistor, SI 

R(s) example-system input command signal (frequency domain) 
r(t) example -system input command signal (time domain) 

s Laplace operator, 1/sec 

T sine -wave period, sec 

t time, sec 

X(s) position (frequency domain), cm 

X position (time domain), cm 

3 

V volume, cm 

p 

|3 bulk modulus, N/cm 

^ damping ratio 

2 / 4 

p hydraulic fluid density, N'sec /cm 

T time constant, sec 

w frequency, rad/sec 

I I zero -to -peak amplitude of sinusoid 

Subscripts: 

A peak-power point 

a actuator 

c coil 

cc command compensator 

cl closed loop 

com command 

discharge 


d 


e error 

eff effective 

full full stroke 

h hydraulic 

max maximum 

n numbers in sequence 

nl saturating nonlinear device 

nh hydraulic natural frequency 

ns servovalve -spool natural frequency 

ol open loop 

p plant 

pa power amplifier 

peak peak amplitude 

pre precompensated command 

r rated 

s spool 

sat saturation 

set set point 

su supply 

t total 

V airflow valve 

0 null 

1, 2. . . numerical subscripts used as noted in text 
Superscripts: 

( )* optimum 

( ) first derivative with respect to time 

( ) second derivative with respect to time 

( ) high-pass-filtered signal 
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APPENDIX B 


ITERATIVE PHASE DETERMINATION FROM LIMIT LINES 

Once the limit lines for a particular system design have been determined, the phase 
response of that system can be approximated by combinir^ the limit-line equations with 
the system linear dynamics. Basically, the system is treated as a linear system with 
the saturating device replaced by a linearized gain In the linear operating region, 

is set to 1. When the system is operating in saturation, the effective gain of the 
saturating device is reduced < 1). 

The actual value of is determined by allowing the system output amplitude to 
equal the limit -line amplitude at that frequency and solvii^ the linear transfer function 
for Since this will always result in a below 1. 0, the solution is implemented 

by an iterative process. That is, is initialized at 1. 0 and is decreased in incre- 
ments until the linear-transfer -function magnitude equals the most restrictive limit-line 
magnitude at the selected frequency. The resulting linear transfer function is assumed 
to be correct, and the phase angle of the system at that particular gain and frequency is 
assumed to be the system phase response. This procedure results in a system gain that 
is a function of frequency and input excitation amplitude . 

The closed-loop linear transfer function for the example system of figure 1 is 

C(s) ^ ^nl^p^e (gl) 

as^ + bs^ + s + 

The computer listing that was used to determine the phase plots of figure 7(a) is shown 
here. 


C ITERATIVE PHASE TECHNIQUE 

C PROGRAM #1 - ASYMPTOTIC LIMIT LINES FOR EXAMPLE SYSTEM 

REAL KE,KP,KNL,MAGC 
PI=3. 1415926 
WN=2.*PI*100. 

WI=2.*PI*5. 

DELTW=WI 

1 WRITE(6,23) 

READ{5,2)KE,KP,FSAT,R,ZETA 

WRITE(7,21) 

2 FORMAT (5F8. 4) 

W=WI 

3 D=W/WN 
KNL=1. 

CMAG=R*KNL*KP*KE/SQRT( (KNL*KP*KE-2 .*ZETA*D*W)**2+W*W*(1 . -D*D)**2) 
IF(W-WN)4,6,6 
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4 CLL=(4./PI)*KP*FSAT/W 

GO TO 8 

6 CLL=(4./PI)*KP*FSAT/(W*D*D) 

8 IF(CMAG-CLL)10, 10,12 

10 GO TO 14 

12 MAGC=R*KNL*KP*KE/SQRT( (KNL*KP*KE-2.*ZETA*D*W)**2+W*W*(1 .-D*D)**2) 

IF(ABS(MAGC-CLL).LT..01)G0 TO 14 
KNL=KNL-.0001 
GO TO 12 

14 THETA=- (ATAN2(W*(1 . -D*D) , (KNL*KP*KE-2. *ZETA*D*W) ) ) 

THETAD=THETA*360 . / (2 . *PI ) 

F=W/(2.*PI) 

WRITE(7,16)F,THETAD,CMAG,CLL,KNL 

W=W+DELTW 

16 F0RMAT(1X,6F12.6) 

IF(F-300.)3,3,18 
18 WRITE(6,20) 

READ(5,22) I 

20 F0RMAT(1X,’ RUN ANOTHER CASE? 1=YES,0=N0 ') 

21 F0RMAT(1X,4X, 'FREQ' ,7X, 'THETA' ,9X, 'MAG' ,9X, 'CLL' ,10X, 'KNL' ) 

22 FORMAT(Il) 

IF(I)24,24,1 

23 FORMATdX,' INPUT KE,KP,FSAT,R,ZETA ') 

24 STOP 
END 


The iterative phase technique is not quite as straightforward when used with a light- 
ly damped system containing a velocity feedback loop (figs. 6 and 7(b)). When the sys- 
tem operates linearly, the velocity feedback increases the damping of the forward-loop 
transfer function. But when the system operates in saturation, this feedback is opened 
during the saturation portion of the cycle. 

The portion of the cycle when the system is in saturation can be found by solving for 
the time at which the error-signal sine wave reaches saturation or 

KjE| = (B2) 


Solving for tj gives 


t. = 1 sin-1 ( sat 

CO 




(B3) 


If the saturation is assumed to be well established, or Fsat/^l^e I E I ) ^ 0* 5, the approxi- 
mation 


sin-1 


vKeEl/ K^E 


(B4) 


26 



is used to simplify equation (B3), giving 


^sat 

Ct)Kg|E 


(B5) 


In terms of the sine-wave period, w can be replaced by 27 t/T. Since tj^ is the time 
from t = 0 to the point where the first satiiration occurs, there will be four of these 
time periods in a cycle. The ratio of the time in linear performance to the period of 
the sine wave is 


4t, 


2F 


sat 


"■Ke E 


(B6) 


To estimate the damping ratio for a system with a velocity loop that is closed only a 
portion of the time, we assume a simple straight-line relationship for the damping ratio. 

?aff = ?ol + «cl - (B’* 


The proportionality constant K is related to the degree of saturation and varies from 
1. 0 for no saturation to 0 for saturation over the entire period. Equation (B6) defines 
the percentage of time that the system is not in saturation and is used to give an approx- 
imate value for K. Therefore, 

(BS) 


This equation provides a good approximation for the effective damping ratio result- 
ing from the velocity feedback loop with the system in partial saturation. However, 
without an exact time-solution at all frequencies, it is difficult to determine the value of 
|E I . It is approximated by a value of | R|/2, which was foimd to give good agreement 
when used to predict the phase of the lightly damped example-system case. Therefore, 
the relationship used to determine the effective damping ratio of the fojrward-loop trans- 
fer function for partial saturation is 


<e£l = ?ol + «cl - 5oI> 


4F 


sat 


UrKe R 


(B9) 


This damping ratio is switched in when the linear response crosses the limit line before 
the saturation gain is decremented. A listing for the program used to generate the data 
of figure 7(b) is presented here. 
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C ITERATIVE PHASE TECHNIQUE 

C PROGRAM #2 - COMPLEX LIMIT LINE FOR EXAMPLE SYSTEM 

REAL KE,KP,KNL,MAGC 
PI=3. 1415926 

1 WRITE{6,23) 

READ ( 5 , 2 ) WN , W I , DELTW , KE , KP , FSAT , R , ZETOL 
WRITE(7,21) 

2 FORMAT (8F8. 4) 

DK=.l 

W=WI 

3 D=W/WN 
KNL=1. 

CMAG=R*KNL*KP*KE/SQRT((KNL*KP*KE-2.*.7*D*W)**2+W*W*(1.-D*D)**2) 
DLL= KP*FSAT/ ( W*SQRT ( D*D*D*D+4 . *D*D* . 04-2 . *D*D+ 1 . ) ) 

CLL=DLL*4./PI 
IF(CMAG-CLL)5,5,10 
5 ZETA=.7 

GO TO 15 

7 ZETA=ZETA+.01 

GO TO 12 

10 ZETA=ZETA0L+(ZETA-ZETA0L)*(FSAT*4. )/(KE*R*PI) 

12 MAGC=R*KNL*KP*KE/SQRT((KNL*KP*KE-2.*ZETA*D*W)**2+W*W*(1.-D*D)**2) 
IF(ABS(MAGC-CLL).LT..01)G0 TO 15 

IF(MAGC-CLL)13,15,14 

13 IF(KNL.EQ.1.)G0 TO 7 
KNL=KNL+DK 
DK=DK/2 

GO TO 12 

14 KNL=KNL-DK 
GO TO 12 

15 THETA=- (ATAN2(W*(1 . -D*D) , (KNL*KP*KE-2. *ZETA*D*W) ) ) 

THETAD=THETA*360 . / (2 . *PI ) 

F=W/(2.*PI) 

WRITE(7,16)F,THETAD,CMAG,CLL,KNL,MAGC,ZETA 

W=W+DELTW 

16 F0RMAT(1X,7F12.6) 

IF(F-300.)3,3,18 

18 WRITE(6,20) 

READ (5, 22) I 

20 F0RMAT(1X,' RUN ANOTHER CASE? 1=YES,0=N0 ') 

21 FORMAT ( IX, 4X, 'FREQ' ,7X, 'THETA' ,8X, 'CMAG' ,8X, 'CLL' ,9X, 'KNL' ,9X, 'MAG- 
IC ,10X, 'ZETA' ) 

22 FORMAT(Il) 

IF(I)24,24,1 

23 FORMATdX,' INPUT WN,WI, DELTW, KE.KP, FSAT, R.ZETAOL ') 

24 STOP 
END 
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APPENDIX C 


DERIVATION OF PEAK-POWER EQUATIONS 

The peak-power level delivered to a hydraulic actuator piston was determined in 
reference 9. The derivation of the peak-power ellipse and the resulting theoretical peak 
values of flow and pressure are repeated here for completeness. 

The output, or load acceleration, can be repressed in terms of the pressure across 
the piston. Substituting = |X^|sin cot into equation (15) and rearranging result in 

A 

-5. P (t) = - |X 1 0)2 sin o)t (C 1) 

M„ ^ ^ 


where 


Pa(t) = Pi(t) - P2^t) (^2) 

When equation (19) is used, the output velocity can be expressed in terms of the servo- 
valve output flow as follows: 


q«(t) 


= X 


X^ I CO cos cot 


(C3) 


Rearranging equations (Cl) and (C3) results in the following two equations: 

— = cos 0)t 

CO 



(C4) 




-sin cot 


(C5) 


Therefore, for sinusoidal motion of the output x^(t), equations (C4) and (C5) show 
that the servovalve output spool flow q„(t) and pressure p_(t) must also be sinusoidal 
but with a specific phase relation. Figure 16 is a typical plot of the phase relation of the 
variables x„(t), q„(t), and p„(t). The valve flow q„(t) leads the output by n/2 radians, 
and the flow and pressure are out of phase by tt/2 radians. When the flow and pressure 
for one complete sine wave are crossplotted on a servovalve characteristic, a load locus 
is generated. An expression for this load locus is now derived. 
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Squaring both sides of equations (C4) and (C5) and adding the results yield 


a2 2 


= 1 





(C6) 


Equation (C6) describes an ellipse and represents the locus that would appear on a 
servovalve ou^ut characteristic plot of flow versus pressure for sinusoidal motion of 
the actuator output position x^(t). 

The largest possible ellipse must be bounded by the maximum spool displacement 
±jX j and the maximum piston pressures ±|P„ | • To get the most out of the 

system, the ellipse described by equation (C6) must lie within this region but be tangent 
to the ±|X I curves at some point. The choice of this tangent point is somewhat 

arbitrary but will be selected as the point where the servovalve transfers maximum 
power to the load. This point is designated as point A in figure 17. The derivation of 
this point is considered in reference 12 and, therefore, only the results are presented 

here. The pressure p. and flow q_ at this peak -power point are as follows: 

3. s 


^a, A o ^su 




(C7) 

(C8) 


where Q^, is defined as the servovalve rated flow that normally appears in the specifi- 
cations. Rated flow is conventionally defined as the maximum flow available from the 
valve with (2/3)Pg^ across the load. 

The derivation for the load locus or ellipse that is tangent to the |X | curve at 

“ max 

point A of figure 17 is now presented. The expression for the servovalve pressure-flow 

characteristic at maximum spool displacement |X | is as follows: 

® max 


where 


% 


= Xo 
s ' s 


max 



(C9) 



(CIO) 
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The generalized expression for the pressure and flow output requirements during sinu- 
soidal output motion is given by equation (C6). After rearranging, this equation be- 
comes 


q 


2 

s 



(Cll) 


To determine a specific expression of the load ellipse that is tangent to the servovalve 
characteristic at point A (fig. 17), we determined and equated the derivatives dq^/dp^ 
of equations (C9) and (Cll). The result is 


CsFsl p / A^ 

max / a 


2VPsu - Pa 

Substituting qg = C |X | - p and rearranging give 

max » 


(C12) 


2(P - P ) 

' su ^a' 



(C13) 


Substituting in the coordinates of point A of figure 17 ((2/3)P_„, Q_) yields an expression 

SU 1 

for the frequency of oscillation for the output locus that passes through point A: 


a>A = 


- Po 
3 ^ 

^vQr 


(C14) 


Substituting equation (C14) into equation (Cll) at the coordinates of point A yields the 
peak displacement |X^ for the locus through point A: 



3v/2q2m^ 


(C15) 


Substituting equations (C14) and (C15) into the general expression of equation (Cll) yields 
the expression for the maximum output (load) pressure and flow locus: 
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(C16) 



Rearranging yields the following: 



(C17) 


This is the expression for the load locus or ellipse that is tangent to the |X | curve 

max 

at point A in figure 17. In terms of supply pressure P_„ and rated flow Q„, its inter- 
cepts with p„ and axes are as follows: 

a. S 


*a, peak g ^ ^su 

(C18) 

^s, peak ~ Qp 

(C19) 


These are the peak piston pressure and flow that the servovalve must supply to enable 
the output load to generate the locus of equation (C17) when the system is operating in a 
sinusoidal manner. The values of peak pressure and flow defined by equations (C18) and 
(C19) then determine the output acceleration and velocity limit criteria, respectively. 
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TABLE I. - EXAMPLE-SYSTEM PARAMETERS 


p 

Example-system transfer -function coefficient, a, sec 2. 533X10” 

Example-system transfer -function coefficient, b, sec: 

At damping ratio ? of 0. 7 2.228X10”^ 

At damping ratio ^ of 0.2 6.366X10"^ 

Force level at saturation, V 8 

Error gain, 10, 5, 1 

Plant gain, 50, 100, 500 




TABLE n. - AIRFLOW -VALVE SERVOSYSTEM PARAMETERS 


Controller: 

Precompensator lead frequency, rad/sec 

Precompensator lead damping ratio, ^ 

Precompensator lag frequency, cug, rad/sec 

Precompensator lag damping ratio, C 

High-pass-fUter time constant, sec 

Integrator low-frequency cutoff time constant, t^, sec 

Error gain, K^, V/V 

Proportional to spool velocity gain, K„_, V/V 

ps 

Integral of spool velocity gain, Kj^, sec“i 

Power amplifier and servovalve coil: 

Power amplifier resistors, 

• 

«2 

R4 



^8 



Power amplifier voltage -to -current gain, V/A . 

Coil resistance, 

Coil inductance, L^, H 

Coil force -to -current gain, K^, N/A 

Servovalve spool dynamics: 

Spool mass, M^, N * sec^/cm 

Spring rate, N/cm 

sp 

Spool spring-mass system damping, N* sec/cm . 
Servovalve orifices: 

Null flow -pressure coefficient, K^p q, cm^N’Sec . . 

Spool diameter, dg, cm 

Orifice discharge coefficient, C^ 

Leakage orifice area at null, cm2 

Hydraulic fluid: 

2 

Hydraulic supply pressure, P^„, N/cm 

2/4 

Density, p, N • sec /cm 

2 

Bulk modulus, /3, N/cm 


. 2ir(315) 
. . 0. 50 

. 2jt(950) 
. . 0. 50 

- . 0. 02 
. . 0.01 
. . 1. 70 

. . 5.04 

1.313X10^ 


. . 0.28 
. 40 200 

. 11 000 
... 127 
. . 2000 
. . 2000 
. 25X10^ 

... 3.0 
1. 0X10"^ 
. . 13.34 

1. 124x10'^ 

. . . 4404 

. . . 0.078 

7. 15X10"^ 
. . . 1.316 
. . . . 0.6 
1. 548X10"^ 

. . . 2068 
8. 175X10"^ 
6. 895X10^ 


Actuator piston: 

Full stroke (zero to peak), cm 0.318 

Area, A_, cm^ 3.200 

» o 

Total hydraulic fluid volume in actuator, V^, cm 28. 68 

Mass of piston and load-valve moving element, M^, N • sec^/cm ... 1. 122X10"^ 

Damping factor, D^^, N* sec/cm 1.016 

Transducer gains: 

Load-valve position transducer, K^, V/cm 15. 75 

Servovalve -spool velocity transducer, V • sec/cm 0. 0472 

Saturation levels: 

Coil-current saturation, A 6.2 

Coil-voltage saturation, eJ^ V 35 



TABLE m. - SUMMARY OF AIRFLOW -VALVE 


LIMIT -LINE EQUATIONS 


Limit line 

Equation 

Value 

Acceleration limit 

jX^lo.2 = 

^ 3My 

5. 561X10® cm/sec^ 

Coil-current limit - 
low frequency^ 


2. 407x10^ cm/sec 

Coil-current limit - 
high frequency^ 


9. 430X10® cm/sec® 

Coil-voltage limit - 
at low frequency^ 

ix |a.=JWe?^ 

^ JrK^jR„ + RjA. 

sp^ c l a 

4. 142X10^ cm/sec 

Coil-voltage limit - 
at high frequency** 

|X |o,^ = ^^^cV,sat 

5. 324X10^ cm/sec^ 


Vhere = irdgC^ y = 3. 221X10^ cm^/sec 





Figure Z. - Block diagram of unity-gain feedback system with saturation in forward loop. 
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Plant gain, 



Figure 2. - Dynamic performance limit lines of example system for various plant gains. 
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amplitude, 

|R| 

9 



Analog mode! 

Linear-nonlinear bound 

Limit line 


J L^ l -I I I I I I J ■ I I I . I ■ 

10 100 
Frequency, Hz 

(c) Plant gain, Kp, 500; error gain, Kg, 1. 
Figure 3. - Concluded. 










System output amplitude, |C|, 





Phase, C(s)/R(s), deg 



Figure 7. - Analog-computer-model phase response and approximate phase analysis. Plant gain, Kp, 100; error 
gain. Kg, 5. ^ 


























Phase angle between load-valve position and command, deg 



Frequency, Hz 

Figure 15. - Phase-angle response of airflow valve compared with iterative 
analysis technique for predicting phase. 



0 tt/2 tt 3tt/2 2tt 

Angular measure, cat, rad 


Figure 16. - Typical position, flow, and pressure 
response during sinusoidal motion. 



1. Report No. 

NASA TP-1488 


2. Government Accession No. 


3. Recipient's Catalog No. 


4. Title and Subtitle 

PREDICTING DYNAMIC PERFORMANCE LIMITS FOR 
SERVOSYSTEMS WITH SATURATING NONLINEARITIES 


5. Report Date 
July 1979 


6. Performing Organisation Code 


7. Author(s) 

John A. Webb, Jr., and Richard A. Blech 


8. Performing Organization Report No. 

E-9903 


9, Performing Organization Name and Address 

National Aeronautics and Space Administration 
Lewis Research Center 
C le veland, Ohio 44135 


10. Work Unit No. 

505-05 


11. Contract or Grant No. 


12. Sponsoring Agency Name and Address 

National Aeronautics and Space Administration 
Washington, D.C. 20546 


13. Type of Report and Period Covered 

Technical Paper 


14. Sponsoring Agency Code 


15. Supplementary Notes 


16. Abstract 

This report demonstrates a technique for predicting the dynamic performance limits of a system 
with saturations. A generalized treatment for a system with a single saturating nonlinearity is 
presented and compared with frequency-response plots obtained from an analog model of the sys- 
tem. Once the amplitude dynamics are predicted with the limit lines, an iterative technique is 
employed to determine the system phase response. The saturation limit-line technique is then 
used in conjunction with velocity and acceleration limits to predict the performance of an electro- 
hydraulic servosystem containing a single-stage servovalve. Good agreement was obtained be- 
tween predicted performance and experimental data. 


17. Key Words (Suggested by Author(s) ) 

Electrohydraulic servosystem; Servosystem; 
Servovalve; Fast-response system; Nonlinear; 
Saturation; Analysis; Control 


18. Distribution Statement 

Unclassified - unlimited 
STAR Category 08 


19. Security Classif. (of this report) 
Unclassified 


20. Security Classif. (of this page) 

Unclassified 


21. No. of Pages 

50 


22. Price* 

A03 


For sale by the National Technical Information Service. Springfield. Virginia 22161 


NASA-Langley, 1979 




National Aeronautics and 
Space Administration 

Washington, D,C. 

20546 

Official Business 

Penalty for Private Use, $300 


V. ' 

fVIASA 


THIRD-CLASS BULK RATE 


\ : 


Postage and Fees Paid 
National Aeronautics and 
Space Administration 
NASA-451 



5 VI 10, A, 072079 S009030S 

dipt of the AIB EOBCE 

AF BEAPONS lABOBATOBl 

ATTH: TECHNICAL IIBBABY (SOI) 

KIBTLAND AFB NM 87117 




POSTMASTER: 


If Undeliverable (Section 158 
Postal Manual) Do Not Return 







